Fourier mode dynamics for the nonlinear Schrodinger equation in one-dimensional 

bounded domains 



O 

(N 



> 

O 

rn 
O 



J. G. Caputoi, N. K. Efremidis^, Chao Hang^.^ 
^Laboratoire de Mathematiques, INS A de Rouen, 
B. P. 8, 76801 Saint-Etienne du Rouvray, France 
E-mail: caputo@insa-rouen.fr 
^Department of applied Mathematics, 
University of Crete, 71409 Heraklion, Crete, Creece 
E-mail: nefrem@tem.uoc.gr 
"^Centra de Fisica Teorica e Computacional, 
Faculdade de Ciencias, Universidade de Lisboa, 
Avenida Professor Cama Pinto 2, 
Lisboa 1649-003, Portugal 
E-mail: chang@cii.fc.ul.pt 
'''Department of Physics and State Key Laboratory of Precision Spectroscopy, 
East China Normal University, 
Shanghai 200062, China 
(Dated: January 18, 2013) 

We analyze the ID focusing nonlinear Schrodinger equation in a finite interval with homogeneous 
Dirichlet or Neumann boundary conditions. There are two main dynamics, the collapse which is 
very fast and a slow cascade of Fourier modes. For the cubic nonlinearity the calculations show no 
long term energy exchange between Fourier modes as opposed to higher nonlinearities. This slow 
dynamics is explained by fairly simple amplitude equations for the resonant Fourier modes. Their 
solutions are well behaved so filtering high frequencies prevents collapse. Finally these equations 
elucidate the unique role of the zero mode for the Neumann boundary conditions. 



I. INTRODUCTION 



The nonlinear Schrodinger equation in 2 or 3 spatial 
dimensions with a focusing nonlinearity has been studied 
intensively because of the collapse phenomenon, see for 
example the reviews of Berge [l| and Sulem Most 
of these studies have been done for an infinite domain, 
for which the equation is invariant by a scaling transfor- 
mation. This symmetry is important to determine the 
conditions for collapse. In many applications like for ex- 
ample a laser propagating in an optical fiber, the domain 
is finite so the boundaries play an important role. A pi- 
oneering study was conducted by Fibich and Merle Q 
for the 2D cubic nonlinear Schrodinger equation. They 
showed that circularly symmetric ground state waveguide 
solutions are stable and that the critical power condition 
for collapse is sharp unlike for an infinite domain. For 
small amplitudes the ground states reduce to the Bessel 
linear modes of the Laplacian. An interesting problem 
is then how these modes exchange energy as the solu- 
tion evolves. This issue is important since some of these 
modes can be filtered out. More generally this energy ex- 
change between linear modes is related to the old problem 
studied by Fermi-Pasta and Ulam (see the first section 
of [4] for subsequent developments). For a chain of os- 
cillators with cubic nonlinearities, a medium amplitude 
Fourier mode initial condition gives rise to a cascade of 
higher Fourier modes and energy fiows back into the ini- 
tial mode. For the cubic nonlinear Schrodinger equation 
in ID we expect a similar recurrence because the underly- 
ing long wave reduction model remains the Korteweg-de 



Vries equation. This recurrence is different from what 
happens in turbulence where there is a one way fiow of 
energy in wave numbers. For example Muraki [5] stud- 
ied this one-way cascade for the Burgers equation. Using 
the Cole-Hopf transformation, he was able to quantify 
the phenomenon. Note finally the study by Zakharov 
et al of wave turbulence carried out on a ID model in 
Fourier space. 

Following a similar approach, in this article we have an- 
alyzed the ID nonlinear Schrodinger equation with cubic 
and quintic nonlinearities on a finite interval with Dirich- 
let or Neumann boundary conditions. We have chosen to 
work in ID to benefit from the Fourier machinery which 
enables to solve the problem fairly easily. A relatively 
small number of Fourier modes are necessary to describe 
well the solution when it is not singular. Another advan- 
tage is that the analysis can be done easier than for the 
2D case. We have studied the stationary solutions and 
their stability. We have also obtained simple models for 
the resonant transfer of energy between Fourier modes for 
the cubic and quintic nonlinearities. As expected there 
is no transfer for the cubic nonlinearity reflecting the in- 
tegrability of the equation. For a quintic nonlinearity a 
resonant transfer exists. The solutions of these reduced 
models compare well to the numerical solutions of the 
partial differential equation. This method of analysis 
of the resonant transfers of energy can be extended to 
higher dimensions. The principle remains the same but 
of course the machinery will be much more complicated. 
From another point of view our study presents the time 
evolution of a solution prior to collapse. Finally we show 
that filtering prevents collapse. 
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The article is organized as follows. In section II we review 
the conservation laws and the Virial identity, the main 
theoretical tool, for both the Dirichlet and the Neumann 
boundary conditions. Section III presents the ground 
states and their stability. Numerical results are shown in 
section IV and explained using a new model of resonant 
energy transfer in section V. We conclude in section VI. 



II. CONSERVED QUANTITIES AND VIRIAL 
RELATIONS 



We consider the one dimensional nonlinear Schrodinger 
(NLS) equation 

«V't + V'xx + |^|'''^ = o (1) 

on a smooth, bounded domain [0, tt] with the homoge- 
neous Dirichlet boundary condition 



= 0) = xjj{x — tt) — 0, 
or the Neumann boundary condition 

i>x{x = 0) = ijjxix = tt) = 0. 



(2) 



(3) 



Here c? is a positive integer with d = 1, 2, and 3 corre- 
sponding to the cubic, quintic, and septic nonlinearities. 

Eq. ([l} with both Dirichlet and Neumann boundary 
conditions admits the following conserved quantities, the 
norm 



(4) 







which is the total power in optics and the Hamiltonian 
The momentum 



d+1 



2d+2 



dx. 



(5) 



n = ^ / (^V: - r^x)dx, 
/o 



(6) 



which is conserved for the infinite domain now has a fiux 

n, = -4[|V'.p]J, (7) 

for Dirichlet boundary condition and 



2d 
d+l 



1^1 



2d+2 



(8) 



for Neumann boundary condition. 

We can analyze the evolution of following integral 
quantities related with the model ([IJ 



h{t) = n^i^x^dx, 

Jq 

hit) - rfi^.p- 



2{d+l 



2d+2 



dx 



d+1 



2d+2 



dx. 



(9a) 



(9b) 



Here, Ji is the variance, which is a common tool for pre- 
dicting collapse of NLS equation solutions in infinite do- 
main and we assume that Ij {j = 1,2) are initially well 
defined. 

For Dirichlet boundary, some algebra leads to 



dt^ 



Jo d+l 



(10) 

which shows that < is a sufficient condition for col- 
lapse. Notice that the right hand side of pH)) can be 
negative even though H > 0. In the infinite domain the 
first term is absent. Therefore, compared with the in- 
finite line, Dirichlet boundaries focus the solution and 
enhance the collapse. 

For the Neumann boundary conditions, we have 



dt^ 



27r 



2d 



2d+2 



8H -4: 



d-2 

rf+T 



IV'I 



2d+2 



dx. 



(11) 



From Eq. ([TT|) we notice that H < is not a suffi- 
cient condition for the collapse. Therefore the Neumann 
boundary can be either reflecting or absorbing and en- 
hance or suppress the collapse. 



III. BOUND STATES 

The time periodic solutions of Eq. ([TJ can be searched 
in the form 



ipiz, x) = u{x) exp{iEt), 



(12) 



where u{x) is a real function and E is the propagation 
constant which is also real. The resulting equation then 
reads 



-Eu + u" + ^ 



with the Dirichlet boundary conditions 

u{x = 0) = u{x = tt) = 
and the Neumann boundary conditions 

Ux{x = 0) = Ux{x = tt) = 0. 



(13) 
(14) 
(15) 



A conserved quantity of Eq. (jl3p can be found by 
multiplying with u'{x) and integrating over x 



2^ ' 2 



1 



2d+2 



2d + 2 

which can be further integrated to give 



dx — ±— I dz 
2 



z\AK + 2Ez 



^d+l 



(16) 

-1/2 

(17) 
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FIG. 1: Phase portrait associated with Eq. 1 II6II for E = —1 (left 
panel) and E = 1 (right panel) for d = 2. The levels presented for 
K am K = 0, 0.1, 1, 10, 20, 30, 40, 50, 60. 



by defining z = u^. The above integral can be expressed 
in terms of Jacobi elliptic functions in both the cases of 
cubic and quintic nonlinearities. However, in the latter 
case the resulting expressions are rather complicated. 

The phase portrait associated with Eq. ([TS]) is shown 
in Fig.[T]for E — ~1 (left panel) and E = 1 (right panel). 
The points at m = (ux — 0) correspond to the solutions 
on the boundaries for Dirichlet (Neumann) case. 

If u <C 1, the nonlinearity does not play an important 
role and the solutions are close to the linear limit of Eq. 
dm, i.e. 



Um{x) = Asin(ma;), 
for Dirichlet boundary condition and 
Um{x) — Acos(ma;), 



(18) 



(19) 



for Neumann boundary condition. Here m denotes the 
index of the wave number m = 0, 1, 2, • • • and A denotes 
the amplitude with A 1. The propagation constant, 
K , and the power are then given by 



En 



K - —42 



P^^A\ 
2 



(20) 



When |u| starts to increase, the propagation constant, 
as well as the form of the solutions, is slightly modified 
due to the effect of nonlinearity, i.e. E — + 6E with 
SE <^ Em- The relationship between the correction of 
the propagation constant 5E and amplitude A (and thus 
P) can be computed by 



5E^ 



(2d + 2)! 



22<i+l[(d+l)!]2 



(2d + 2)! 



^<i2'i+i[(d+l)!]2 



(21) 



When M ~ 1, the nonlinearity has a strong effect and we 
have to resort to numerical methods to find the solution. 
We also notice that Eq. (IT3l) with the Neumann bound- 
ary condition admits constant solution, i.e. u = Aq with 
E = Al'^> 0. 

In Fig. [21 we solve Eq. ([T^ by using a shooting method 
to construct stationary modes for both Dirichlet and 
Neumann boundary conditions and find authentic law for 
P{E). The stationary solutions u{x) for d = 2 are shown. 
The initial conditions are taken by the solutions ((T5)) and 
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FIG. 2: The stationary solutions u{x) for d = 2. The initial 
conditions are taken as sin(r?i2:) and cos(rria;) naturally satisfying 
the Dirichlet and Neumann boundary conditions, respectively, with 
m = 1 [the upper two panels in (a)] and m = 2 [the lower two panels 
in (a)]. The dependence of P on _E are shown in (b) and (c) for 
m = 1 and m = 2, respectively. The total power P\e^oo ~^ 2.755 
for m = 1 and P|b_»oo ~^ 5.447 for m = 2. The red line denotes 
the value of Per by II25I I. 



(IT^ . We see that the family of solutions with m = 1 bi- 
furcates from zero at E = Ei = — 1 while the family of so- 
lutions with TO = 2 bifurcates from zero at E = E2 — —4. 
Close to the bifurcation points the P — E curves follow 
the relation (|2ip . P{E) is monotonically increasing with 
the growth of E while dP{E) / dE\E^oo 0. 

The linear stability of stationary solutions can be nu- 
merically computed by solving the linearized eigenvalue 
problem. Particularly, we assume that 



iEt 



(22) 



where a, j3 are small perturbations which are propor- 
tional to exp(— iAt). The coupled eigenvalue problem 
then reads 



LiB = -i\A, L2A = i\B, 



(23) 



where Li = —dxx + E — u{x) and L2 = —dxx + E — 
{2d+l)u{x)'^'^ . The growth rate is defined as max[Im(A)] 
at which an unstable solution will grow. We find that the 
family of solutions with to = 1 is always stable against 
linear perturbations whereas only a very narrow stability 
region close to the bifurcation point exists for the family 
of solutions with to = 2, 3, 4 • • • . 

In the case of infinite domain, Eq. (IT^ admits the 
localized soliton solution 



u=[{d+ l)£;]i/2''sechi/'' \d^/E (x - I 



(24) 



where the solution is centered in the center of the domain 
X = 7t/2. The total power for d = 2 is a constant 



p — 

^ r.r — 



^dx = ^^12 = 2.72, (25) 
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i.e. it is independent on E. The solutions with power 
P < Per disperse during the propagation, whereas if 
P > Per the solutions collapse. A similar analysis can 
be carried out for the case of d = 3 leading to 



(a) 



Pr. 



\u\'^dx 



25/37ri/^r(7/6) 1 
3i/2r(2/3) £;i/6' 



(26) 



which is expressed in terms of F functions and depends on 
E. Although the soliton solutions do not satisfy the NLS 
equation on a bounded domain with specific boundary 
conditions, they are particularly useful as limiting cases 
of solutions. 

The solution (|24l) can be considered as good approx- 
imations for the solutions satisfying both Dirichlet and 
Neumann boundary conditions in the limit E ^ oo. This 
is because in this latter limit the maximum intensity in- 
creases as E^^'^ while the pulse width decreases as E^^^. 
The narrowing of the pulse makes the soliton tails, as well 
as their derivatives, are almost zero on the boundaries. 
The critical power for collapse Pd does not change with 
the growth of E. 



IV. THE NUMERICAL SIMULATIONS 

To understand the dynamics of the Fourier modes, we 
relied heavily on numerical solutions of the NLS equa- 
tion on the interval [0,7r]. The equation was solved us- 
ing the split-step Fourier method where the linear part 
is advanced using the sine or cosine Fourier transforms 
respectively for Dirichlet and Neumann boundary condi- 
tions. The details of the numerical implementation are 
given in the Appendix A. 

We first consider Dirichlet boundary condition. The 
solution can be expanded in sines as 



E 



(t) sin mx. 



(27) 



For d = 1, there is no collapse for Eq. ([T]). As expected 
the evolution of a sine initial condition gives rise to a 
cascade of modes. For ci(0) = 2 and Cj^i(O) = we 
observe a cascade to C3 and C5 with maximum amplitudes 
max(c3) — 0.7 and max(c5) = 0.15 with the other modes 
being insignificant. For 03(6) = 2 and Cj^3{0) = we get 
almost no cascade. 

Now let us compare the outcomes for d — 1 with 
the initial conditions 'ip{0,x) =^ sin(a:;) + 2sin(3x) and 
"0(0, x) = sin(x) -|- sin(3a;). The time evolution of the 
mode amplitudes are shown in Fig. [3] The mode ampli- 
tudes fluctuate in a fairly narrow range around an average 
value. This range decreases even more for smaller initial 
amplitudes as shown in the right panel. We will explain 
these effects in the next section. 

Let us now turn to the quintic case (c? — 2) with 
Dirichlet boundary condition. In Fig. 21 we show the 
time evolution of solution (jlSp with m = 1 and differ- 
ent values of amplitude A. No collapse is observed when 




FIG. 3: Time evolution for a cubic nonlinearity and Dirichlet 
boundary conditions of the sine mode amplitudes Cm{t), m = 1, 3, 5 
starting from two different initial conditions ^{0, x) = sin(a::) -|- 
2sin(3x) (left panel) and ip{0,x) = sin(a:) -|- sin(3a:) (right panel). 



A (P) is small whereas collapse occurs when A > 1.3 
{P > 2.65). A collapsing solution is shown on the last 
row of Fig.glfor A = 1.31. The right panels of Fig.|3]show 
the Fourier spectra. We notice that only the odd modes 
are excited. Actually, we can explain that the cascade 
of Fourier modes for Eq. ([T]) starting with a particular 
mode g (g = 0, 1, 2, • • • ) is always restricted to the modes 
q(2n — 1) (n = 1, 2, 3, • • • ) irrespective of the boundary 
conditions. In other words, we can expand the solutions 
of Eq. (dl) as 



ip{t,x) 



00 



C,(2„-l) 



i[g(2n - l)a;] 



(28) 



for the initial condition ■0(0, x) — sm{qx). For the Neu- 
mann boundary conditions the sines should be changed 
to cosines. Details are given in Appendix B. In Fig. [51 we 
show the recurrence of the solutions and spectrum cas- 
cade at different times when no collapse occurs. We will 
explain this phenomenon in the next section. We also 
notice that for large propagation constant, the solutions 
appears in the form of the hyperbolic secant function. 

In Fig. [51 we show the time evolution of the Fourier 
amplitudes for a initial condition (|18p with m — 1 and 
increasing amplitude A. Each panel calculated for a sin- 
gle amplitude corresponds to different regions of similar 
behaviors for the modes. When < A < 0.5, there is 
only the mode m = 1. The other modes are insignificant. 
When 0.5 < A < 1.0, we observe two modes m = 1 and 
3. When 1.0 < A < 1.2, we observe three modes m — 1, 
3, 5. When 1.2 < A < 1.3, four modes m — 1, 3, 5, and 
7 are observed. The larger the amplitude of the initial 
condition, the more modes are excited. When A > 1.3 
collapse occurs. This is a much faster mechanism than 
the recurrence. The energy travels very suddenly from 
the low frequency modes to the higher frequency modes. 

For the Neumann boundary conditions, the zero mode 
plays a special role. For example, it does not give rise to 
a spectrum cascade. We will see below how it couples to 
the other modes. Apart from this, the simulation results 
for the Neumann boundary conditions are similar to the 



5 



(a) 
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FIG. 4: Left panels: time evolution of the maximum of \ip\ for 
a quintic nonlinearity and Dirichlet boundary conditions for three 
different initial conditions, from the top to the bottom, VCO, a^) = 
l.Osin(x) (P = 1.57), 1/^(0, x) = 1.3sin(x) (P = 2.65) and ■ip{0,x) = 
1.31sin(a;) (P = 2.70). The right panels show the corresponding 
Fourier spectra at time 4 = 4. Notice the collapse occurring at 
t Si 0.7 in the last row. 




(b) 



FIG. 5: Snapshots of the wave profiles and Fourier spectra for a 
quintic nonlinearity and Dirichlet boundary condition. The times 
shown are t = 0, 0.6, 1.2 and 1.8 respectively in panels a, b, 
c and d. The initial condition is ^l>{0) = 1.3sin(a;) (P = 2.65). 
The recurrence of the wave profile and the spectrum cascade is 
evident. The profiles shown in (c) and (g) follow the hyperbolic 
secant function \2A^ . 



ones for the Dirichlet boundary conditions. 



V. EVOLUTION OF THE RESONANT 
FOURIER MODES 

In the interval of existence of solutions, one can use 
the Fourier sine or cosine series depending whether we 
have Dirichlet or Neumann boundary conditions. In both 
cases we have a complete basis for functions satisfying the 
boundary conditions. We will first consider the Dirich- 
let boundary condition and examine cubic and quintic 
nonlinearities. The results that we will obtain for these 
are very similar for the Neumann boundary condition un- 



(a) 


A=1.0,P=1.57 




1 


2 3 A 
t 


(c) 




A=1.2,P=2.26 







(b) 



1 

0.5 
0^ 



A=1.1,P=1.90 




FIG. 6: Time evolution of the Fourier amplitudes for quintic non- 
linearity and Dirichlet boundary condition for different initial con- 
ditions. In panels a, b and c, ci(0) = A = 1., 1.1 and 1.2, respec- 
tively. The panel d shows max(|cm|) for m = 1, 3, 5, 7 and 9 as a 
function of the initial amplitude A. 



less the zero mode is involved. This specific case will be 
addressed at the end of the section. 

For the Dirichlet boundary conditions, we expand ^ in 
a sine Fourier series 



i^i^^t) = X! c„i{t)sm{mx), 



where the are given by 



Cm{t) — — i/j{x,t) am{nx)dx. 

TT 



(29) 



(30) 



We find that if coUapse occurs, it happens due to the 
divergence of the Fourier series. Indeed, the conservation 
of the power implies for any m' 



(31) 



because of Parseval's relation, 



£\^\^dx^^r. 



Therefore the Fourier coefficients are bounded and the 
only source of collapse is the divergence of the series. 

Substituting the expansion ([^ into Eq. ([!]) , we obtain 
the coupled equations of Fourier amplitudes. For the 
quintic nonlinearity {d = 2) we get 



iCq - q Cq + 
= 0, 



CkCiCmc'^c*{klmnp\q) 



k.l,m.n,p 



(32) 



where 



(klmnplq) = / sin(fcx) sin(Za;) sin(mx) sin(na;) 
Jo 

X sin(px) sin{qx)dx. 
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Eqs. (|32)) can be simplified by the transformation Cq 



yielding 



TT 



xe 



akaia,na'^a*{klmnp\q) 

k,l,m,n,p 
i(k^ +f -n'^ ~q^)t _ g 



(33) 



Similarly, the coupled equations of Fourier amplitudes for 
the case of septic nonlinearity (d = 3) read 



— akaiamana*a*a*{klmnpqr\s) 



k,l.m,n.p.q.r 



xe 



0. 



(34) 



Note that in Eqs. (l33l) and (|34p most terms are rapidly 
rotating and average out to zero. Only the ones such 
that fc2 + /2 ^ _ _ p2 _ g2 ^ for (1221) ( fc2 + /2 + 

m^ + n^— p^ — q^ — r^ — — for i.e., the resonant 
terms will contribute to the long term dynamics of aq. 

A detailed study can be carried out of the dynamics of 
the Fourier coefficients for different types of nonlineari- 
ties. We have used the Maple software [1] to identify all 
the resonant terms in the equations for the mode ampli- 
tudes. For the cubic nonlinearity {d = 1) the amplitude 
equation equivalent to (j33p reads 



TT 



/ , akaia„^ 

k,l,m 



{klfn\j) 



0. 



(35) 



Taking into account the resonance condition, Eq. ([35 
turns into the following equations 



laj + cij \V 



where V 
solution 



is conserved. Eq. (|36l) admits the 



a,- e 



(37) 



An obvious implication is that d/dt{\aj\'^) = so that 
there is no transfer of energy from one mode to another. 
This is what we have seen in the numerical results in 
Fig. [21 Over a short time the modes oscillate in a pe- 
riodic fashion, however their average over a long time is 
constant. As expected no collapse will occur in the cubic 
NLS equation irrespective of the initial total power. In 
fact the integrability of the NLS equation on the whole 
line, related to the existence of a Lax pair, has been 
shown to carry over to the case of a finite domain with 
Dirichlet boundary conditions 

The case of the quintic nonlinearity {d — 2) is more 
complicated. For simplicity, we consider a solution con- 
sisting of three modes, i.e., to = 1, 3 and 5, 



4'{t,x) = flie cos(a;) -I- 036 
-I- a5e"*^^*cos(5a;). 



cos(3a;) 



(38) 



Then Eq. ([33l) gives rise to the following coupled resonant 
amplitude equations 



idi + oi 



+ ^(l«3|^ + |a5n 
"9 



8 



alaial = 0, (39a) 



ias + as 



13 



2|„ |2 



asl'V - —lasl" + 3\ai\'\a5 



"9 



a^af as = 0, (39b) 



16 



ios -I- 05 



9 



\a5\^V-^\a,\^ + 3\a,\'\as\' 



+ ^{\ai\' + \asn 



16 



a*^a^ = 0, (39c) 



0, (j = 1, 2, • • • 00) (36) defining /, 



where V = |aip + |a3p + [asp. It is easy to check that 
Eqs. 1211) satisfy the condition ^(|ai p-Hjaap-^lasp) = 0, 
i.e., is a conserved quantity. The last terms in equa- 
tions represent the mixing between the three modes and 
therefore the intensity of each mode is not conserved as 
in the case of d = 1. At this point note that if we had 
included ay in the description we would have had the 
additional resonant terms 

I |2 2 2| |2 * 2| |2 * 

If I ay I << 1 then their contribution would be very small. 
Another point is that if we had included the even modes 
in the description we would have obtained the extra reso- 
nant term a2a4^a5 in the equation for oi. This indicates 
that the modes 2 and 4 couple through mode 5. 

Eqs. (|39p can be written into a more compact form by 
{j — 1, 3, 5). Here, Ij represent the 
intensity of each mode satisfying V = Ii + 13 + 1^- Then, 
we obtain 



i 3 3/2 1/2 . . 



h = 



--/1/3' /g' sm^ 
3/2ri/2,- /J 



(40a) 
(40b) 
(40c) 



where 9 = — 26'i + Wt, — 9^ with 9j — argOj. We notice 
that the driving terms on the right hand side of Eqs. (HOI) 
are proportional to the intensities of the modes and their 
phases. In addition, we have the constraints 



3-^1 + 9^3= Ml, (41a) 

^^3 + \h = /i2, (41b) 
4 8 

-^h - = Ml - M2 = M3, (41c) 

where fXi, (z = 1 — 3) are constants of the motion. Using 
the above relations, the dynamics can be reduced to the 
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equations for Ii and 6. They are 

A = ^/i/^/'/g/'sin^, (42a) 

e = 1(13/2 + 3/f + 7/| + 21/1/3 + 27/1/5 + 12/3/5) 

-'ihll'^h^'\ (42b) 

together with the constraints (HJ). 

In order to check the equations for the resonant Fourier 
modes, we compare the solutions of the reduced equa- 
tions and the NLS. In Fig. [7] we show the time evo- 
kition of the solutions of Eqs. (|40l) and ^ by using 
the Runge-Kutta method and split-step Fourier method, 
respectively. The initial conditions are the same. As 
expected the amplitudes of the Fourier modes for ^ 
present fast periodic oscillations. However over a long 
time interval the solutions of the reduced equations and 
the full partial differential equation match well, support- 
ing the validity of the reduced model. 

At this point, let us consider the Neumann boundary 
conditions for which there is the additional zero mode. 
For the cubic nonlinearity c? = 1, the evolution of the 
Fourier modes follows (|36p so that there is no resonant 
energy exchange between the modes /^ = 0, j = 0, 1, . . . . 
For the quintic nonlinearity d = 2, the situation is more 
interesting. Assuming a solution containing the three 
modes i = 1, 3, 5, we obtain evolution equations identical 
to (|39)) except that the signs of the resonant terms are 
reversed. We then obtain the same final equations 
except that the evolution of the phase 6 is reversed. 

Assuming a solution containing the three first modes 
i = 0,1,2 , we obtain evolution equations of the form 
(p9| but with no terms outside the brackets. This means 
that again, no resonant transfer of energy exists between 
modes. If the third mode is added to the expansion, 
new terms appear outside the brackets. The evolution 
equations are 





h flo [. . . 


1 3 , 2 * 

J + -00010203 : 


= 0, 


(43a) 


idi - 


hai [... 


1 3 2 *2 

J + -0002 03 = 


0, 


(43b) 


10,2 - 


h 02 [. . . 


] -1- ^aoOia203 : 


= 0, 


(43c) 


10.3 - 


h 03 [. . . 


] + l^faial = 


0, 


(43d) 



where the . . . terms in the brackets are all real. Following 
a similar procedure as above, the modal energies Ij — 
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FIG. 7; Comparison between the solution of the reduced model 
and the numerical solution of the NLS for quintic nonlinearity and 
Dirichlet boundary condition. The top left (right) panel shows 
the time evolution Ii{t), I'iit), and l5{t) {d{t)) for the reduced 
model I I42II . The bottom left (right) panel shows the corresponding 
evolutions for the NLS. The initial conditions are taken as /i = 
0.64, I3 = 0.36 and h = 0.16 and 9 = -tt. 



lojp evolve as 



/o 


3 J jl/2 J jl/2 . a 

= 2^0// sm6l. 


(44a) 


A 




(44b) 


/2 


= 3/o/y'/2/3^'sin0. 


(44c) 


is 


3 r rl/2 T rl/2 • a 


(44d) 



where 6 ~ —26*0 + Oi + 262 — 63. As above one can then 
reduce the problem to two equations, one for Iq and one 
for 9. We do not write these equations because they are 
cumbersome. The interesting fact is that one needs the 
4 modes 0-3 present in order to see this resonant transfer 
of energy. If one of them is missing there is no energy 
transfer. This particular feature of Neumann boundary 
conditions changes the route for collapse for the Dirichlet 
case and the Neumann case. 



VI. FILTERING FOURIER MODES PREVENTS 
COLLAPSE 

The projection of the nonlinear Schrodinger equation 
on a a finite number of Fourier modes will yield ampli- 
tude equations that are well behaved and do not exhibit 
collapse. As we have seen collapse is related to a sud- 
den energy flow to high frequencies. Thus, one can ar- 
rest the collapse by filtering the high Fourier modes and 
therefore preventing the sudden energy flow to high fre- 
quencies. Physically this can be done by introducing a 
nonlocal absorption in the model ([T]). This conclusion is 
also available for the models with septic nonlinearity and 
Neumann boundary condition. 
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FIG. 8: Left panels: Time evolution of the maximum of |i/>| for 
quintic nonlinearity with different initial conditions. Right panels: 
The corresponding Fourier spectra at t = 4. The initial conditions 
are taken as ip{0,x) = 1.4sin(x) (P(0) = 3.08) in the first row 
and '0(0,2;) = 1.4sech(x — ^) (P(0) = 3.56) in the second row. 
We compare the results with and without filtering the high Fourier 
modes designated by the black and red lines, respectively. The 
collapse is efficiently prevented by filtering in both cases. 

In Fig. |8] we show the time evolution of solution with 
or without filtering high Fourier modes. The initial con- 
ditions are taken as a sine wave (the first row) and a 
hyperbolic secant pulse (the second row). In both cases, 
the collapse is efiiciently prevented by filtering. Because 
of the absorption of the higher modes, there is a small 
loss (less than 10 percent) of the total power. This can 
be decreased by including more lower frequency modes. 



VII. CONCLUSION 

We have analyzed and solved numerically the one di- 
mensional nonlinear Schrodinger equation on a finite in- 
terval with Dirichlet or Neumann boundary conditions. 
A preliminary analysis reveals that i7 < is sufficient 
for collapse in the Dirichlet case but not for the Neu- 
mann case. The bound states have been computed. The 
first non trivial one corresponding in the linear limit to 
sin a; (resp. coscc) for Dirichlet (resp. Neumann) bound- 
ary conditions is always linearly stable as opposed to the 
higher order modes whose window of stability is very 
small and reduces as the order is increased. 

We have solved the partial differential equation for the 
cubic and quintic nonlinearities. In the cubic case there 
is no resonant transfer of energy between Fourier modes, 
while it is present in the quintic case. Identifying res- 
onant terms in the evolution equations of the Fourier 
modes, we have written reduced systems. Their evolu- 
tion is in excellent agreement with the solutions of the 
NLS, even close to collapse. For the Neumann boundary 
conditions the Goldstone mode plays a particular role as 
it couples the modes m = 1, 2, 3. For the Dirichlet case, 
there is no coupling between the first three modes. Fi- 



nally note that this model reduction can be extended to 
higher dimensions and other systems like a cylindrical 
waveguide. The machinery would be more complicated 
but the overall method remains the same. 
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Appendix A: Numerical procedure for solving the 
ID NLS © 

The ID nonlinear Schrodinger equation ([IJ with 
Dirichlet or Neumann boundary conditions is solved as 
usual by splitting the linear and nonlinear part of the 
operator. The linear part 

iut + u^x = 0, (Al) 

is such that 

where u is the sine or cosine Fourier transform of m, to 
satisfy the boundary conditions. The solution of the non- 
linear stepping is the standard one 

u{2dt) = e*l"('**)l"''**u(dt). 

The numerical implementation is done in Matlab Q 
and the solution is evaluated at discrete points u„,n = 
1, . . . A'^. The sine and cosine Fourier transforms for the 
linear step (IA1|) are then done using the discrete sine and 
cosine Fourier transforms. 

For the Dirichlet boundary conditions, we use the dis- 
crete sine Fourier transform 

u{k) = f2^HM^^), k = l,...N, (A2) 

71—1 

and inverse discrete sine Fourier transform 

u(n) = ^u(fc)sm(-^-^), n^l,...N. (A3) 

k=l 

The Neumann boundary conditions are trickier to im- 
plement because one needs to use the discrete cosine 
Fourier transform 

uik) ^ w{k) f2 nin) cos( ^^^^ ~ ^^^^ ~ ), = 

71—1 

(A4) 
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with 



w(l) = l/ViV, w(fc) 



2 < fc < 



. The inverse discrete cosine Fourier transform is 
u{n) = y w{k)u{k) cos( ^^^" ~ \y ~ n = l,...iV. 



fe=i 



(A5) 

The number of discretisation points was chosen to be 
TV = - 1 or iV = - 1 ^-^^Yi a step dt = 10"*. The 
L2 norm was checked during the computation and it is 
conserved up to 10~i° in absolute value. 



Appendix B: Parity argument for the mode cascade 



In order to prove Eq. ([28l) . let us consider the case of 
d = \, the cases of d = 2, 3 can be proved in a similar 
way. Eq. ^ with d — 1 can be written into the form 



tpit + di, x) = ip{t, x) + i 
+ \ij{t,x)\'^'ip{t,x)] dt. 



ld^ip{t,x) 
2 dx^ 



(Bl) 



Now we assume that the initial condition is taken as 
V'(0,a;) = cos(gx) = i(e*«"= + e-'?^). Substituting the 
initial condition into Eq. (jBip . we have 



1 



1 



3iqx 



3e^«^ + 3e- 



4 

■3iqx ^ 



dt. 



-iqx ^ 

(B2) 



As we see, after a short time interval dt, the solution 
'ip{dt,x) can be expressed as ip{dt,x) = CqCos{qx) + 
C3q cos{3qx). We can repeat this process and obtain that 
^it,x) = Em=(2«-l)gCmC0s(mx) (u = 1,2,---). In 
other words for c? = 1 we immediately obtain modes 1, 
3 and 5 while for d = 2 we immediately obtain modes 1, 
3, 5 and 7, etc. Finally, notice that for Neumann bound- 
ary condition there is no cascade to higher Fourier modes 
starting with the m = mode. 



Appendix C: Evolution of the resonant Fourier 
modes for the Neumann boundary conditions 

For the Neumann case, we expand ^ in a cosine Fourier 
series 

00 

ip{x,t) = ^ Cm{t) cos{mx), (CI) 

m=0 

where the Cm are given by 



1 r 

co{t) = - t}j{x,t)dx, 



(C2a) 



2 r 



Cn(t) = - / i)(x,t)cos{nx)dx (n 7^ 0). (C2b) 



Substituting the expansion (|C1|) into Eq. ([T]), we ob- 
tain the coupled equations of Fourier amplitudes for the 
case of quintic nonlinearity [d = 2) 



ICQ 



— y CkCiCmC^cl{klmnp\0) = 0, (C3a) 

TT ^ — ' 



iCq-q^Cq + - y CkClCjnCnC*{klmnp\q) 



0, {q ^ 0) 



(C3b) 



where 



{klmnp\q) — / cos(fca;) cos(Za;) cos(ma;) cos(nx) 
^0 

X cos(pa;) C08{qx)dx. 
Eqs. (jC3p can be simplified by the transformation Cq — 



yielding 



■q + — akaiamana*{klmnp\q) 

TT ^ — ' 

-/? 2 I .2 I 2 2 2 2\i 

+/ +m -n -p -9 )* — Q 



(C4) 



where Uq — 1 for q — Q and Uq — 2 for q ^ Q. Similarly, 
the coupled equations of Fourier amplitudes for the case 
of septic nonlinearity {d ~ 3) read 

ids H akaiamanO*pa*^a*{klmnpqr\s) 

y_^-ih(k'^+l'^+m^+n^-p^-q^-r^-s'^)t ^ ^(-(5^ 

where Ug = 1 for s = and (jg = 2 for s ^ 0. 

For simplicity, we consider the solution psp . Then, 
Eq. (|C4p turns into the following coupled equations 



icti -I- oi 



9 



^|aipP-H|«^|4+3|,3|2|„^|2 



+ 8(l«3|' + |a5n 
"9 



103 -I- 03 



9 



+ -alalal = 0, (C6a) 
13, 



-\a,\^P--\a,\^ + 3\a,\'\a,\' 



+ g (|ai|^ + |a5 
"9 



-alafa.^O, (C6b) 



13, 



-\a,\'P--\a,\^ + 3\a,\'\a,\' 



a3 



_3 
16 



-^af a| = 0, (C6c) 



where V = \ai\^ + jaap -I- |a5p. These are exactly the 
same as p9|) except that the terms outside the brackets 
have the opposite signs. 

As done above, eqs. (IC6p can be written into a more 
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compact form by defining Ij = \aj\ {j = 1, 3, 5), i.e. together with the constraints (|4T 

h^^h,4^'ll/^sm0, (C7a) 
/3 = ^/i/|/'/5'/'sin0, (C7b) 
h = -lhl!^'ll^^ sine, (C7c) 



where 9 = -2611 + 303 - 65 with = arga^. Eq. ([C7)) 
can be further written into the form 

A = -7/1/3^^/5^^ sin 6*, (C8a) 



4 

1 

' 4 



-(13/2 + 3/2 + 7/2 + 21/1/3 + 27/1/5 + 12/3/5 



3/i/|/'/5-^/'], (C8b) 
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